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Bose-Einstein-condensed gases in external spatially random potentials are considered in the frame of a 
stochastic self-consistent mean-field approach. This method permits the treatment of the system properties 
for the whole range of the interaction strength, from zero to infinity, as well as for arbitrarily strong disorder. 
Besides a condensate and superfluid density, a glassy number density due to a spatially inhomogeneous compo- 
nent of the condensate occurs. For very weak interactions and sufficiently strong disorder, the superfluid fraction 
can become smaller than the condensate fraction, while at relatively strong interactions, the superfluid fraction 
is larger than the condensate fraction for any strength of disorder The condensate and superfluid fractions, 
and the glassy fraction always coexist, being together either nonzero or zero. In the presence of disorder, the 
condensate fraction becomes a nonmonotonic function of the interaction strength, displaying an antidepletion 
effect caused by the competition between the stabilizing role of the atomic interaction and the destabilizing role 
of the disorder With increasing disorder, the condensate and superfluid fractions jump to zero at a critical value 
of the disorder parameter by a first-order phase transition. 

PACS numbers: 03.75.Hh, 03.75.Kk, ()3.75.Nt, 05.30.Jp, 05.70.Ce 



I. INTRODUCTION 

Physics of dilute Bose gases is usually studied (see Refs. (H, 
S m B HI S 01) for asymptotically weak interactions, when 
the Bogolubov approximation |H, |^ is applicable. In the 
case of strong interactions, the most reliable techniques are 
purely numerical, such as Monte Carlo simulations fiol [Til 
Mill [3, [m [IB, [I3l- In the majority of experiments with 
trapped atoms, interactions are rather small 0, 0, S S 0] cor- 
responding to values of the gas parameter much smaller than 
one. However, now it has become possible to vary the interac- 
tion strength in a wide range by using the Feshbach resonance 
techniques (see Refs. iSliSl)- Thus, in experiments with ^^Rb 
atoms il9l |20l (2l|], the value of the gas parameter 0.8 has 
been reached. Large values of the scattering length and, re- 
spectively, strong effective interactions can also be achieved in 
quasi-one-dimensional and quasi-two-dimensional configura- 
tions due to t he g eometric resonance 1,22. ,23:1 (see discussion 
in Refs. US Him HI). 

Recently, based on the idea of representative statistical en- 
sembles 1271! as applied to Bose systems with broken gauge 



symmetry 
veloped |3; 



a self-consistent approach has been de- 
I32i HH] for treating Bose-condensed sys- 



tems with arbitrarily strong interactions. This approach was 
shown ySjl to reproduce the weak-cou pling exp ansions of Bo- 
golubov (§, ^ and Lee-Huang- Yang |l34i l35l |^^, while si- 
multaneously being in good agreement with numerical Monte 
Carlo simulations for strong interactions. 

A fundamental feature of any uniform Bose system with 
arbitrarily strong interactions is the appearance, at low tem- 
peratures, of a Bose-Einstein condensate and, simultaneously, 
of superfluidity. In these systems, the condensate fraction, 
no, and the superfluid fraction, rig, always coexist, both being 
nonzero below the condensation temperature. Though there 
is no simple relation between these fractions, the superfluid 



fraction is always larger, rig > hq. 

When a uniform Bose system is subject to the action of an 
external spatially random field, the relation between the con- 
densate and superfluid fractions could change. Thus, Huang 
and Meng |37] considered a Bose-condensed system in a ran- 
dom external potential. They treated the case of asymptoti- 
cally weak interactions and asymptotically weak disorder As- 
suming that their results could be formally extended to strong 
disorder, they suggested that there can exist the so-called Bose 
glass phase, in the sense that there is a Bose-Einstein conden- 
sate, no 7^ 0, but there is no superfluidity, Ug = 0. Weakly 
interacting Bose gas in the presence of disorder was also the- 
oreticafly studied in Refs. ||38i [39i |40i |4 ij Hj, kj. |44tl and the 
experiments demonstrating the localization of atomic matter 
waves were performed recently ll45ll46ll47ll . 

The arising Bose glass phase, if any, would be of high 
importance for the experiments with ^He-filled porous me- 
dia |48]. A porous material can be mimicked well by an ex- 
ternal random potential. This is because each pore represents 
an external local potential. At the same time, since pores en- 
joy random properties, being of different sizes, shapes, and 
being randomly distributed in space, they do form for a Bose 
system a kind of a spatially random potential |48]. However, 
Monte Carlo simulations i49tl as well as numerical calcula- 
tions in the frame of the random-phase approximation fsd^, 
accomplished for a Bose system with strong disorder, though 
with asymptotically weak interactions, revealed no Bose glass 
phase. But maybe this phase could appear when both disorder 
and interactions were strong? In the present paper we shall 
give our answer to this question. 

In Ref. ifsTll . a self-consistent stochastic mean-field approx- 
imation for Bose systems in external random potentials, al- 
lowing the treatment of arbitrarily strong interactions and ar- 
bitrarily strong disorder, have been developed. It has been 
found that the disordered system contains several particle frac- 
tions whose relative size between and 1 defines the main 
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system properties. There exists the fraction no of condensed 
atoms and there is the fraction 12^ of normal uncondensed 
atoms. There also appears the fraction no of a glassy compo- 
nent. An important role is played by the anomalous average 
(7, whose absolute value \a\ quantifies the relative density of 
pair-correlated atoms lf52ll . Finally, there exists the superfluid 
fraction Ug. All these quantities are defined by the solution of 
a system of nonlinear equations, whose exact analysis would 
require numerical calculations. It is the aim of the present 
paper to give an analysis of the behavior of these various frac- 
tions under varying strengths of interactions and disorder in 
the domain where tiq 7^ 0. Throughout the paper we use the 
system of units, where h = 1 and fc^ = 1- 



II. MAIN DEFINITIONS AND NOTATIONS 

First, we must define the system to be considered. The 
Hamiltonian energy operator has the standard form 



H 



V/(r) 



2m 



-e(r) 



ip{r) dr 



+ -f I ^t(r)V,t(r)V(r)V(r)dr, 



(1) 



in which (r) is the Bose field operator, ^ (r) is an external 
random potential, and the interaction strength is 



$0 



0,8 

4tt — 

TO 



(2) 



with the scattering length and the atomic mass to. The 
random potential is centered around zero, so that its stochastic 
average vanishes, 

((e(r)»=0. (3) 
The stochastic correlation function has the general form 

mmr')))=R{r-r'). (4) 

For what follows, it is important to distinguish between the 
stochastic averaging over the random-field distribution, which 
is denoted through the angular double brackets ((. . .)), as in 
Eqs. (|3]l, (IDi, and the statistical averaging over the quantum 
degrees of freedom, which, for an operator A, is denoted as 



{A)h = Ti-pA , 



(5) 



where p — p[H] is a statistical operator having the Gibbs form 
with a grand Hamiltonian H. The total averaging, given by the 
simple angular brackets (. . .), includes both the stochastic as 
well as the quantum averaging. 



(i) ^ (( )) = ((TrM» . 



(6) 



The appearance of a Bose-Einstein condensate implies the 
spontaneous gauge symmetry breaking, which can be realized 
by the Bogolubov shift |53, 54] of the field operator 



where ?7(r) = (V'(r)) is the condensate wave function and 
ipi{r) is the Bose field operator of uncondensed atoms. The 
field variables 77(r) and V'i(r) are treated as two independent 
variables, orthogonal to each other. 



i]*{r)'4!i{r) dr — , 



(8) 



which excludes the double counting of the degrees of freedom. 
The quantum-number conservation condition 



(9) 



defines ?7(r) as the system order parameter, equal to the aver- 
age (-0(r)). 

These two field variables obey two normalization condi- 
tions. The condensate function is normalized to the number 
of condensed atoms 

No^ J h(r)P dr, (10) 

while the number of uncondensed atoms 

Ni - (iVi) (11) 
normalizes the number operator 

Ni = J^lir)Mr)dr (12) 

for uncondensed atoms. The system stability is guaranteed by 
minimizing the grand thermodynamic potential 



(13) 



under the constraints of the two normalization conditions (fTOl l 
and ( fTTl i. which requires the use of two Lagrange multipliers, 
so that the grand effective Hamiltonian in Eq. (fT3T l is 



(14) 



with H being the energy operator ([T]i under the Bogolubov 
shift dTji. The form of the grand potential (fTSl l corresponds to 
the quenched disorder 

In view of the zero-centered random potential, satisfying 
Eq. (|3]l, the system can be considered as uniform on average, 
such that 



ri{r) = (po = N^/V) 



(15) 



The field operator of uncondensed atoms can be expanded 
over plane waves. 



ik-r 



(16) 



The operators ak in the momentum representation define the 
normal average 



ip{r) '0(r) = r/(r) + ViCr) , 



(7) 



Uk = {a\ak) 



(17) 
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and the anomalous average 



(18) 



The major quantities to be studied are the densities of dif- 
ferent components. The condensate density po = p — pi is 
expressed through the total average density 



p = N/V = po + Pi 



(19) 



where N is the total number of atoms, and through the density 
of uncondensed atoms 



Pi 



The anomalous average 



-E 

V ^ 



V ^ 



CTk 



(20) 



(21) 



gives the density | cri | of pair-correlated atoms . In the 
presence of random fields, there appears an additional impor- 
tant quantity, the density of the glassy component 



PG 



^-^E((i(«'»)^i')) 



(22) 



fc/O 



whose definition is analogous to the Edwards-Andersen order 
parameter for spin glasses fsS*]. 

Superfluidity is characterized by the superfluid fraction 



1 d ,^ 

hm — • P,, 



3mN v^o 



(23) 



defined as the fraction of atoms nontrivially responding to the 
velocity boost with the velocity v, under w = |v| ^ 0. Here 



P + TOVTV 



(p = E^4a/cj (24) 



is the total momentum of the moving system, and the averag- 
ing in Eq. ( |23] ) implies that with the Hamiltonian of the mov- 
ing system 



Hj, = H 



N + Y,{^-^)aUk , (25) 



where 



N = No + Y^ alak ■ 

k^O 

It can be shown (see, e.g., Refs. iH Isill ) that definition ( |23] | 
yields 

2Q 



where 



n„ = 1 



Q 



3T ' 



2mN 



(26) 



(27) 



is the dissipated heat per atom. 



III. STOCHASTIC MEAN-FIELD APPROXIMATION 

Since our aim is to treat arbitrarily strong interactions and 
disorder, we cannot neglect any part of the total grand Hamil- 
tonian ( fT4] i. For example, if we would omit the terms of the 
third and fourth order, containing the products of three and 
four operators or aj^, as well as the third-order term, in- 
cluding the product aj,ap£^k-p, we would come to the Bo- 
golubov approximation, used by Huang and Meng [37], and 
many others, which allows the consideration of only asymp- 
totically weak interactions and disorder Contrary to this, we 
shall retain all terms of the Hamiltonian, using the stochastic 
mean-field approximation of Ref. ifsill . This approximation 
was previously shown to give a very accurate description of 
different statistical systems with stochastic fields, as is sum- 
marized in Refs. ifSa, IstI Issll . The stochastic mean-field ap- 
proximation for Bose systems with random fields has been 
described in full detail in the recent paper ifsil] . Therefore in 
the present work, we limit ourselves by mentioning only the 
principal points of this approximation and by reviewing the 
resulting formulas that are necessary for the following anal- 
ysis. In the following we shall choose the condensate wave 
function r/ as real without restriction of generality. 

The third- and fourth-order terms of the Hamiltonian with 
respect to the products of the operators and a[, are treated 
by means of the Hartree-Fock-Bogoluboy approximation, 
similarly to the case without random fields i30ll3lii32il33ll . A 
special care is taken with regard to the third-order term con- 
taining the random field ^k- To this end, we use the fact that 
there are two types of averaging, the stochastic averaging, as 
in Eqs. ([3]), (|4]i, and (l50l l. and the quantum statistical averag- 
ing, as defined in Eq. Q. Let us introduce the random quan- 
tity 



ak 



[akjH 



(28) 



which is the quantum average (lU of the operator Ofc. The 
random variable (l28T l is not 0, though its stochastic average 

(K)) = (afc) = 

is zero because of condition (|9]l. Then, using notation ( |28] l. we 
accomplish a mean-field-type decoupling for the third-order 
term 



(29) 



where only the quantum averaging is involved, but no stochas- 
tic averaging is taken. Keeping here the stochastic averaging 
unapproximated makes it possible to consider any strength of 
disorder 

The following important step is the use of the nonuniform 
and nonlinear, with respect to the random variable ^j., canon- 
ical transformation 



flfe = UkOk + V-kb-k 



LUk + mc^ 



(30) 
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containing a new random variable (pk to be defined later Here, 



^fc + gfc 
2m 



L^k - £k 



2£t 



mc 



2\2 



(31) 



The latter equation, for ej., can be represented as the Bogol- 
ubov spectrum 



1.2 



(32) 



however with the sound velocity c differing from that of the 
Bogolubov form. Rather, c is given here as the solution to the 
equation 



(33) 



The random variable ipk, introduced in the canonical trans- 
formation (l30l l. has to be chosen so that to simplify the total 
Hamiltonian. If the variable ipk satisfies the Fredholm equa- 
tion 



_ VNq 

'^'^ V V ^ uj. 

p#0 



i.k-p^p 



then the Hamiltonian acquires the simple form 

H = EB+y] ekblbk + ^/% (fo 



(34) 



(35) 



k^O 



in which the quantum variables bk and 5^ are separated from 
the random variable ipk, and where the first term is a c-number 
quantity 

k^O 

The last term in Hamiltonian ( |35] ) is obtained by using 
Eq. ( [34l l, with the expression for tfo defined by the equation 



fo 



V 



mc 



2 ■ 



The relation between the random variable ak, defined in 
Eq. ( |28] ), and the random variable ipk, satisfying Eq. ( |34] |. fol- 
lows from Eq. ( l30l l, from where 



(36) 



With Hamiltonian ( 1351 1, it is straightforward to calculate the 
normal average ([TtI i, which gives 



nk 



^ coth f 
2efc \2TJ 



((KP» 



(37) 



and the anomalous average ( fTSl ), yielding 



i^coth(|l)+((Kn). (38) 



The last terms in Eqs. dJTl l and (l38T l are caused by the random 
potential. According to relation ( |36l ), we have 



((KP» 



(39) 



The density of uncondensed atoms ( [20] l consists of two terms. 
Pi = PAf + Pg , (40) 
the first of which is the density of normal uncondensed atoms 



1 



sk y^Tj 



1 



(2^ 



(41) 



and the second term is the density ( |22] ) of the glassy compo- 
nent 



PG 



{{\fk\')) 



dk 



(uJk + rnc2)2 (27r)3 
The anomalous average (ISTT i also is a sum 

CTl = (Tat + PG 

of the term 



tTjV 



coth 



dk 



(2^) (27r)3 



(42) 



(43) 



(44) 



and of the glassy density ( |42] ). 

The superfluid fraction ( |26] l is expressed through the dissi- 
pated heat dZTl i. For the latter, we find 



g = g 



N 



where the first term 



dk 



8mp J sm\\{ek/2T) (27r)3 



(45) 



(46) 



is the heat dissipated by normal uncondensed atoms and the 
second term 



1 



2mp J ek{ijJk + mc^) 



coth 



dk 



\2T) (27r)3 



(47) 



is the heat dissipated by the glassy component caused by the 
random potential. 



IV. 5-CORRELATED RANDOM POTENTIAL 

To proceed further in practical calculations, we must spec- 
ify the type of random potential. For this purpose, we take the 
Gaussian (5-correlated random potential with the local corre- 
lation function 



i?(r) - i?o ^(r) . 
Then, by means of the Fourier transformation 



C(r) 



(48) 



(49) 
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the stochastic correlator (|4|l reduces to 

mQ) = hpRov 



(50) 



where V is the system volume. The calculation of the stochas- 
tic average ((\(fik\'^)), involving the method of self-similar fac- 
tor approximants ll59[l60ll . can be done as has been thoroughly 
explained in Ref. jSlll . 

For what follows, it is convenient to deal with dimension- 
less quantities. We introduce the notation for the condensate 
fraction 



no ^ po/p , 
the normal fraction of uncondensed atoms 

nN = Pn I P , 

the glassy fraction 

= PgIp , 

and for the dimensionless anomalous average 

(T = crpf/p ■ 



(51) 



(52) 



(53) 



(54) 



The interaction strength is chai^acterized by the gas param- 
eter 



The dimensionless temperature is 



(55) 



(56) 



The strength of disorder is quantified by the disorder parame- 
ter 



~ Anp^/^ 

Finally, we define the dimensionless sound velocity 

s = racj p^^^ . 



(57) 



(58) 



We shall consider the case no ^ 0, where the gauge sym- 
metry is spontaneously broken. The case of the unbroken 
gauge symmetry with no = will be considered in future 
work. The sound velosity s is then well-defined and given by 
the solution to Eq. ( [33T l, which in the dimensionless quantities, 
and remembering expression ^ for the interaction strength, 
takes the form 



4717(1 — nAT + cr) 



(59) 



The condensate fraction ( BTI i can be found from Eqs. ( fT9] l 
and ( |40] |. which can be reduced to the equation 



riQ ^ ni<s + uq — \ 



(60) 



For the normal fraction ( |52] ). taking into account Eq. ( |4TI ). we 
have, for no 7^ 0, 



37r2 



coth 



3 



2t 



/•oo 



1/2 



1 



dx 



(61) 



The anomalous average (|54l), with Eq. ( l44l i and the dimen- 
sional regularization ||5l|5i|], becomes 



2s2 



7r3/2 
X coth 



7no 



2V27r2 



1)1/2 



2t 



dx 



For the glassy fraction (l53T l. we obtain from Eq. (|42]) 



no = 



v{\ - un) 



(62) 



(63) 



Finally, for the superfluid fraction ( |26] l, employing Eqs. ( |27] l, 
gill, dUll, and gTli, we find 



(64) 



eV^nHJo VTTx^ ai-nh^{s2x/2t) ' 

In the following we shall concentrate on the ground-state 
properties of the system, corresponding to the limit of zero 
temperature t — > 0. In this case, the Bose system without 
disorder would be completely superfluid, rig — 1, but the con- 
densate fraction is depleted by interactions. External random 
fields deplete both the condensate and superfluid fractions, so 
that no < 1 and n^ < 1. What would be the behavior of these 
fractions when varying the interaction and disorder strengths, 
that is, the gas parameter (ISSl l and the disorder parameter jSTli? 

At zero temperature, the normal fraction ( |6TI ) reduces to 



UN 



while the anomalous fraction 



3^ ' 



2s' 
7.3/2 



becomes 



(65) 



(66) 



with nc and no being defined from Eq. (I63I I and normaliza- 
tion d60] l. respectively. The superfluid fraction (|64] i at zero 
temperature takes the form 



1 



no 



(67) 



Equations ( |59] ). ( |60] |, ( |65] ). ( |66] |. and ( l67b determine their so- 
lutions as functions of two variables, the gas parameter ( l55b 
and the disorder parameter ( |57| i. 

Let us, first, consider the asymptotic behavior of the so- 
lutions. If the disorder parameter v is finite and 7 tends to 



6 



zero, there are no physical solutions, which corresponds to the 
stochastic instability of the ideal Bose-condensed gas |51]. 

For a finite disorder parameter and very strong interactions, 
such that 7 oo, we have 



1 /tt 



64 V 9 



5\ 1/3 



with the limit 
Then the normal fraction (|65]| becomes 



-2\l/3 



riAT ~ 1 



TT 

64 



1 



with the anomalous fraction ( |66] l being 



47 \7'^ 
For the glassy fraction (|63] |. we find 



Normalization (|60] | gives the condensate fraction 



no 



6473 



Equation Wh yields the superfluid fraction 



y3 ' 

(68) 
(69) 

(70) 

(71) 
(72) 

(73) 
(74) 



As we see, strong interactions tend to destroy a Bose-Einstein 
condensate, suppressing no, but increase the superfluid frac- 
tion ris. In this limit, Ug ^ uq. 

If the gas parameter is kept finite, but the disorder strength 
is getting asymptotically weak, such that — > 0, then the 
sound velocity is 

s ~ So(l - bu) , (75) 

where the limit 

So = (37r2)i/3a (76) 
depends on the interaction strength through the equation 



47 



and the value of b is given by the equation 

b 



(2 + a3)yTr^_,M,3J ^(1_,3) 



(77) 



So • 
(78) 
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FIG. 1: Condensate fraction no (solid line), superfluid fraction Us 
(dotted line), and glassy fraction no (dashed line) as functions of the 
dimensionless disorder parameter v for very weak interactions, with 
the gas parameter 7 — 10~^. 
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FIG. 2: The same fractions, no, n^, and no, as in Figure [T] but for 
slightly larger interactions, with the gas parameter 7 = 0.001. Here 
the superfluid fraction becomes smaller than the condensate fraction, 
when disorder increases, but Us is never smaller than the glassy frac- 
tion no- 



The solutions to Eqs. ( fTTl ) and ( fTsT l for weak interactions 
(7 0) are 



(97r)i/6 



y/2 + 



16 



(97r)2/3 



b y^7 + 0(7^/^) 



(79) 
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FIG. 5: Fractions no, Us, and no at a relatively large gas parameter 
7 = 1. The condensate and glassy fractions are strongly suppressed. 



FIG. 3: Fractions no, n^, and no for the gas parameter 7 — 0.1. The 
superfluid fraction becomes lower than the condensate fraction at the 
disorder parameter ly ~ 0.9. 
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FIG. 4: Fractions no, Us, and no for the gas parameter 7 — 0.5. The 
superfluid fraction is always higher than the condensate fraction. 
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FIG. 6: Normal fraction um (solid line) and anomalous fraction a 
(dashed line) as functions of the disorder parameter f for weak inter- 
actions, with the gas parameter 7 = 0.1. The anomalous fraction a 
is about three times larger than the normal fraction njv. 



while the anomalous fraction becomes 



and for strong interactions (7 00), Eqs. ( |77] ) and ( fTSl ) yield 

TT 



1 



19273 
1 



1344(97r)i/6 \j3 



1 



(80) 



With the interaction strength varying from zero to infinity, the 
value of a increases from to 1, so that < a < 1. But the 
value of b remains small for all interactions, 6^1. For finite 
7, but asymptotically weak disorder, when 1/ ^ 0, the normal 
fraction is 



UN ~ a^{l — 3biy) , 



(81) 



2sl 

7r3/2 



1 - q3 + (28 - 49a3)6so ^ 



(82) 

where a and b are the same as above. For the glassy fraction, 
we have 



no 



1 - a3 
7so 



The condensate fraction behaves as 

1 - a3 - 21a36so 



no — 1 — a 



7so 



(83) 



(84) 
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FIG. 7: Normal, um, and anomalous, cr, fractions as functions of the 
disorder parameter u for the intermediate strength, with 7 = 0.5. 
The values of a and um are close to each other, though, a > un- 
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FIG. 8: Normal, njv, and anomalous, a, fractions for rather strong 
interactions, with the gas parameter 7=1. Here, contrary to Figs.|6] 
and [7] the normal fraction becomes larger than the anomalous one, 
though they are close to each other. 



FIG. 9: Dimensionless sound velocity s as a function of the disorder 
parameter 1/ for different gas parameters. Each line is marked by the 
corresponding 7. 




0.05 0.1 0.15 0.2 0.25 



FIG. 10: Condensate fraction (a) and superfluid fraction (b) as func- 
tions of the gas parameter 7 for =0 (1), 0.25 (2), 0.5 (3), 0.75 (4), 
1 (5), 1.25 (6). 



For the superfluid fraction, we obtain 

4(1 -a3) 

ric ~ 1 V 

21so 



(85) 



These formulas, to Jeading order in 7 and v, coinside with 
the results of Ref. |37]. In the absence of disorder, the whole 
system would be superfluid, but never completely condensed 
for any finite interactions. Consequently, in the limit of ^ 
0, we have Us > uq. 

With increasing disorder, there occurs a first-order phase 
transition at a value h'c — i^dl), when the system discontin- 
uously transforms to a phase with unbroken gauge symmetry. 



At the point i^c, the fractions no, Ug, no, and a jump to 0, 
while fiN jumps to 1. The behavior of the fractions no, Ug, 
and na as functions of the disorder parameter ly for some se- 
lected increasing values of the gas parameter 7 are shown in 
Figures [T]|5] 

The case of very weak interactions, with 7 = 10^^, is illus- 
trated in Fig.[T] According to Eqs. ( [84l i and ( [85] ), we know that 
at asymptotically weak disorder, when — > 0, the conden- 
sate fraction is smaller than the superfluid fraction, no < Ug. 
But Fig. [T] shows that with increasing disorder the superfluid 
fraction becomes smaller than the condensate fraction. This 
occurs at a rather small value i^, which is too close to to be 
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FIG. 11: Phase diagram in the (7, !/)-plane. The solid line is the 
line of the first-order quantum phase transition between the phases of 
broken and unbroken gauge symmetry. The dashed line indicates the 
values of v and 7 at which the condensate fraction takes the maximal 
value (see Fig. 1 10b). 



noticeable in the figure. Also, we see that disorder suppresses 
the superfluid fraction so that it becomes not merely smaller 
than no, but eventually even smaller than the glassy fraction 
no- 

Increasing the interaction as we go through Figs. |2]|5] 
strengthens superfluidity. Figure |2] for 7 = 0.001, demon- 
strates that, even though there still exists a small value v, 
where the inequality no < Us changes for Us < no, the su- 
perfluid fraction remains now always larger than the glassy 
fraction, n^ > n^. 

When increasing interactions further, say, to the gas param- 
eter 7 = 0.1, as in Fig. [S] the point, where the inequality 
Us > TiQ changes for Ug < no, moves to the right, getting 
closer to the phase transition point i/c- The glassy fraction no 
is always lower than both no and n^. 

In Fig.m for 7 = 0.5, the superfluid fraction is now always 
larger than the condensate fraction, which distinguishes this 
figure from the three previous ones. The condensate fraction 
is yet substantially larger than the glassy fraction. 

Figure |5] emphasizes how strong interactions, with 7=1, 
favor superfluidity, while suppressing both the condensate 
fraction and the glassy fraction. The latter two fractions be- 
come rather small, but the superfluid fraction is close to one. 
This also shows that the system can be practically completely 
superfluid, having a tiny condensate fraction, as it happens in 
liquid ^He. Thus, in Fig. |5] the condensate fraction is about 
5%, though the superfluid fraction is almost 100%. 

The anomalous fraction a and the normal fraction n^r for 
relatively weak interactions, with the gas parameter 7 — 0.1, 
are plotted in Fig.|6l As is seen there, cr is about 3 times larger 
than UN, which stresses the fact that a cannot be neglected. 

For the intermediate interaction strength, the normal and 
anomalous fractions are close to each other, as is shown in 
Fig.|2]for 7 = 0.5. The anomalous fraction is yet larger than 
the normal one. 

When interactions become rather strong, as in Fig. [8] for 
7 = 1, then the normal fraction surpasses the anomalous one. 
But, anyway, is yet close to cr. Varying the disorder pa- 
rameter does not have much influence on the values of un and 
a. 

The dimensionless sound velocity s as a function of the dis- 



order parameter 1/ for different interaction strengths is illus- 
trated in Fig.|9l As it should be, the larger the gas parameter, 
the larger is the sound velocity. Stronger interactions stabi- 
lize the system, increasing the critical value of the first- 
order transition. The sound velocity slightly diminishes with 
increasing disorder 

The condensate and superfluid fractions as functions of 
the interaction strength for different disorder parameters are 
shown in Fig. (TO] As it has been emphasized earlier ifsill . 
the ideal uniform Bose-condensed gas is stochastically un- 
stable. Finite interactions stabilize the system against weak 
disorder But increasing disorder makes the system unstable, 
when the latter transfers through a first-order phase transition 
to a phase with unbroken gauge symmetry. The jumps of the 
condensate and superfluid fractions in Fig. [TOl correspond to 
the phase transition. The superfluid fraction increases mono- 
tonically with the increasing interaction strength. But the re- 
markable fact is that, for nonvanishing disorder parameter, the 
condensate fraction is not a monotonic function of the inter- 
action strength. With increasing gas parameter 7, the conden- 
sate fraction first increases, reaches the maximum, and then 
decreases. Thus, there exists the effect of antidepletion, when 
the increasing interactions result in the rise of the condensate 
fraction. This effect is due to the presence of disorder, which 
tends to destabilize the system, while the interaction stabilizes 
it. The competition between the two tendencies leads to the 
nonmonotonic behavior of no, which is seen in Fig. [TOb . The 
line of the maxima of no in the (7, i/)-pIane is presented as a 
dashed line in the phase diagram in Fig. [TT| 

The line i^c(7) of the first-order phase transitions is drawn 
in Fig.[TT] At the point 7 = 0, corresponding to the ideal Bose 
gas, the phase transition is of second order. However, the ideal 
Bose-condensed gas is stochastically unstable, and for any in- 
finitesimally small disorder parameter v it is destroyed, un- 
dergoing the phase transition to the state with no = 0. Below 
the line Vc{l), there is the superfluid phase, with uq 7^ 0, 
Us ^ 0, nc ^ 0, (T 7^ 0, and nAr < 1. Above this 
line, one has the phase of unbroken gauge symmetry, where 
?io = Us = no = (7 = 0, while nAr = 1. The limit of Vc{'^), 
for 7 tending to infinity, is (Stt^)^/'^. The phase transition 
caused by the increasing disorder is an example of a quantum 
phase transition. 



V. DISCUSSION 

A detailed analysis of the properties of a Bose-condensed 
system at zero temperature in an external random potential 
has been presented. The disorder potential is modelled by 
the Gaussian uncorrected disorder. The strength of disor- 
der as well as the strength of interactions can be arbitrary. 
The system contains several fractions of particles, the conden- 
sate fraction no, superfluid fraction ns, normal fraction n^v, 
anomalous fraction ct, and the fraction of a glassy component 
uq. The behavior of these fractions as functions of the gas 
parameter and the disorder parameter was investigated. The 
ideal Bose-condensed gas is stochastically unstable, since any 
infinitesimally weak disorder destroys it, transferring it to the 
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normal state. Finite interactions stabilize the system. Increas- 
ing disorder leads to a first-order phase transition between 
the superfluid and normal phases. At asymptotically weak 
disorder, such that v ^ 0, the superfluid fraction is always 
larger than the condensate fraction, Ug > riQ. But increasing 
disorder, under very weak interactions, can invert the latter 
inequality, when the superfluid fraction becomes lower than 
the condensate fraction, Ug < uq. This is in agreement with 
the Monte Carlos simulations |49], where it was noticed that 
sufficiently strong disorder can suppress the superfluid frac- 
tion making it smaller than the condensate fraction, provided 
that interaction strengths are very weak. However, at suffi- 
ciently strong interactions, the superfluid fraction gets larger 
than the condensate fraction for all disorder parameters below 
the phase transition point i^c- To our knowledge, no numerical 
simulations have been accomplished, when both the disorder 
as well as interaction strengths would be strong. 

The superfluid and condensate fractions were found always 
to coexist. It may occur that no > Ug or no < Ug, but they are 
nonzero or zero simultaneously. There is no pure Bose glass 
phase, when Ug would be zero, while no is nonzero, though 
the glassy fraction no, induced by disorder, is always present. 

Although disorder suppresses superfluidity, Ug never be- 
comes exactly zero, as one might conclude from the calcu- 
lations for asymptotically weak disorder fs?!] . The pure Bose 
glass phase does not occur in the considered model. 

The unusual effect of antidepletion was found, when in- 
creasing interactions can increase the condensate fraction in 
the presence of disorder This effect is caused by the compe- 
tition between the disorder destabilizing the system and the 
interactions, which stabilize the latter. As a result, in the pres- 
ence of disorder, the condensate fraction becomes a nonmono- 
tonic function of the interaction strength. 

The change in the behavior of the condensed fraction no, 
normal fraction n^v, glassy fraction no, and the superfluid 
fraction n^ results from a competition between the interaction 
potential and the external random potential. These two causes 
act on the fractions in a different way. The increasing inter- 
action always depletes the condensate, but increases the su- 
perfluid fraction. By depleting the condensate, the interaction 
increases the normal fraction n^r. At the same time, strength- 
ening disorder increases the glassy fraction and depletes the 
condensate. The competition of all of these, sometimes con- 
tradictory, governs the overall behaviour of the fractions. 

The origin of the phase transition, occurring under the in- 
creasing disorder, can be understood as follows. Recall that 
the ideal Bose-condensed gas is absolutely unstable with re- 
spect to any infinitesimally weak random perturbations ifsTll . 
Finite interactions do stabilize the Bose-condensed gas. But 
this stabilization can survive only until a finite strength of dis- 
order, when again the system loses stability and transforms 
to the state where the gauge symmetry is not broken. The 
point is that disorder destroys coherence that is pertinent to 
Bose-Einstein condensation. By destroying coherence, disor- 
der moves the system to a state with no Bose condensate. 

The calculations in this paper have been done only for zero 
temperature. This is because, first, it has been necessary to 
understand the behaviour of the system under two varying pa- 



rameters, the interaction strength 7 and the strength of disor- 
der V. Including temperature makes the problem dependent 
on three parameters. This would essentially complicate the 
consideration making it necessary to resort to mainly numer- 
ical calculations. We plan to present the details of these cal- 
culations in our future work. But for low temperatures, the 
obtained results still do hold. Including temperature just leads 
to more condensate depletion and enhancement of the normal 
fraction. 

It is worth emphasizing that when analyzing the behaviour 
of the fractions, we always keep in mind the normalization 
condition (60), according to which the condensate fraction no, 
normal fraction n^v, and the glassy fraction no are added to 
1. However, the explicit relation between the condensate and 
superfluid fractions is not known, because of which the latter 
does not enter any simple normalization condition, except that 
< < 1. 

In conclusion, it is important to discuss the possibility of ex- 
perimental observation of the effects described in the present 
paper. Standard experiments are accomplished with trapped 
atoms. The inclusion of a trapping potential in our theory 
would complicate numerical investigation. However, there are 
two cases, when the results of our consideration could be di- 
rectly applicable to experiments. First, the homogeneous pic- 
ture provides a reasonable approximation for wide traps, and, 
second, it gives a good description of the situation at the center 
of a trap, even if the trap edges are rather sharp. This becomes 
possible because of the known fact that the local-density ap- 
proximation allows for a quite accurate description of trapped 
atoms |1, 2, 3, 4], and the uniform case serves as a starting 
point for the local-density approximation. 

Keeping in mind the local-density approximation, when 
close to the trap center the system is almost uniform, we 
must deal with the gas of atoms with the positive scattering 
length, since a homogeneous gas with attraction is known to 
be unstable iH 12. HI 01 ■ In experiment, one can also realize 
Bose-Einstein condensation of atoms with negative scattering 
length, provided that the atoms are trapped and their number 
does not exceed the critical value Nc- A simple formula for 
the critical number Nc, giving rather accurate estimates for 
harmonic traps can be represented 16111 as 

jY _ Ixlyh 

^" V2 |a,|(/2+Z2+;2) ' 

where la is the oscillator length in the a-direction and the ag 
scattering length. A trapped atomic cloud, with a negative 
scattering length, can be stable only when N < Nc- This 
case requires a separate investigation. In the present paper, 
we have considered a large system with the number of atoms 
not bounded from above. This is why we have assumed from 
the beginning that the scattering length is positive. 

The value of the scattering length can be varied in a wide 
range, for instance, by means of the Feshbach resonance tech- 
niques 13, UMl- It would be interesting to check in exper- 
iment the behavior of the system in a fixed random poten- 
tial, when the interaction strength is varied. Such a situation 
would correspond to Fig. (TO] When diminishing the scatter- 
ing length, that is, diminishing the gas parameter 7, we would 
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come to the boundary of stability of the system. Recall that, 
in the absence of interactions, the Bose-condensed system is 
stochastically unstable, such that any weak random potential 
destroys the condensate. This phenomenon of stochastic in- 
stability was analyzed in detail in Ref. |51]. In order to un- 
derstand, why this phenomenon occurs, it is sufficient to re- 
member that the ideal uniform Bose-condensed gas is unsta- 
ble even in the absence of any random potential, which can 
be easily demonstrated by calculating the system compress- 
ibility and finding out that the latter diverges in the absence of 
interactions ll4l l28ll . 

Finally, the random potential of the type similar to that con- 
sidered in the present paper can be created in experiment, 
e.g., by employing the optical speckle techniques ll45[ 14^ l47h . 
These techniques allow for an efficient regulation of the prop- 
erties of the formed random potential. It is possible to orga- 
nize a frozen random distribution, independent from the time 
variable. It is also feasible to regulate the correlation length 



characterizing the spatial properties of the speckle random- 
ness. When the correlation length is much smaller than the 
healing length, the effective random potential can be repre- 
sented as being J-correlated, which has been assumed in the 
present paper. At the same time, we recall that the general 
theory of Ref. ifsill is applicable to random potentials with 
arbitrary correlation length, although for finite-length correla- 
tions, calculations would be essentially more complicated. In 
this way, it looks quite feasible to check the predictions of the 
suggested approach in experiments with atomic Bose gases 
confined in wide traps. 



Acknowledgments 

This work has been supported by the SFB/TR 12 "Symme- 
tries and universality in mesoscopic systems". 



[1] C. J. Pethick and H. Smith, Bose-Einstein Condensation in Di- 
lute Gases (Cambridge University, Cambridge, 2002). 

[2] L. Pitaevskii and S. Stringari, Bose-Einstein Condensation 
(Clarendon, Oxford, 2003). 

[3] P. W. Courteille, V. S. Bagnato, and V. I. Yukalov, Laser Phys. 
11, 659 (2001). 

[4] V. I. Yukalov, Laser Phys. Lett. 1, 435 (2004). 

[5] J. O. Andersen, Rev. Mod. Phys. 76, 599 (2004). 

[6] K. Bongs and K. Sengstock, Rep. Prog. Phys. 67, 907 (2004). 

[7] V. I. Yukalov and M. D. Girardeau, Laser Phys. Lett. 2, 375 
(2005). 

[8] N. N. Bogolubov, J. Phys. (Moscow) 11, 23 (1947). 
[9] N. N. Bogolubov, Moscow Univ. Phys. Bull. 7, 43 (1947). 
[10] M. H. Kalos, D. Levesque, and L. Verlet, Phys. Rev. A 9, 2178 
(1974). 

[11] D. M. Ceperley, Rev. Mod. Phys. 67, 279 (1995). 
[12] S. Giorgini, J. Boronat, and J. Casulleras, Phys. Rev. A 60, 5129 
(1999). 

[13] J. L. DuBois and H. R. Clyde, Phys. Rev. A 63, 023602 (2001). 
[14] C. C. Moustakidis and S. E. Massen, Phys. Rev. A 65, 063613 
(2002). 

[15] J. L. DuBois and H. R. Clyde, Phys. Rev. A 68, 033602 (2003). 
[16] N. Boninsegni, N. V. Prokofiev, and B. V. Svistunov, 

Phys. Rev. E 74, 036701 (2006). 
[17] S. Pilati, K. Sakkos, J. Boronat, J. Casulleras, and S. Giorgini, 

Phys. Rev. A 74, 043621 (2006). 
[18] R. A. Duine and H. T. C. Stoof, Phys. Rep. 396, 115 (2004). 
[19] S. L. Cornish, N. R. Claussen, J. L. Roberts, E. A. Cornell, and 

C. E. Wieman, Phys. Rev. Lett. 85, 1795 (2000). 
[20] J. L. Roberts, J. P Burke, N. R. Claussen, S. L. Cornish, 

E. A. Donley, and C. E. Wieman, Phys. Rev. A 64, 024702 

(2001). 

[21] N. R. Claussen, E. A. Donley, S. T. Thompson, and C. E. Wie- 
man, Phys. Rev. Lett. 89, 010401 (2002). 

[22] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 

[23] D. S. Petrov and G. V. Shlyapnikov, Phys. Rev. A 64, 012706 
(2001). 

[24] E. B. Kolomeisky, T. J. Newman, J. P. Straley, and X. Qi, 

Phys. Rev. Lett. 85, 1146 (2000). 
[25] M. D. Lee, S. A. Morgan, M. J. Davis, and K. Burnett, 



Phys. Rev. A 65, 043617 (2002). 
[26] G. S. Astrakharchik, D. Blume, S. Giorgini, and B. E. Granger, 

Phys. Rev. Lett. 92, 030402 (2004). 
[27] V. I. Yukalov, Phys. Rep. 208, 395 (1991). 
[28] V. I. Yukalov, Phys. Rev. E 72, 066119 (2005). 
[29] V. I. Yukalov, Int. J. Mod. Phys. B 21, 69 (2007). 
[30] V. I. Yukalov, Phys. Lett. A 359, 712 (2006). 
[31] V. I. Yukalov, Laser Phys. Lett. 3, 406 (2006). 
[32] V. I. Yukalov and H. Kleinert, Phys. Rev. A 73, 063612 (2006). 
[33] V. I. Yukalov and E. P Yukalova, Phys. Rev. A 74, 063623 

(2006) . 

[34] T. D. Lee and C. N. Yang, Phys. Rev. 105, 1119 (1957). 
[35] T. D. Lee, K. Huang, and C. N. Yang, Phys. Rev. 106, 1135 
(1957). 

[36] T. D. Lee and C. N. Yang, Phys. Rev. 112, 1419 (1958). 
[37] K. Huang and H. R Meng, Phys. Rev. Lett. 69, 644 (1992). 
[38] A. V. Lopatin and V. M. Vinokur, Phys. Rev. Lett. 88, 235503 
(2002). 

[39] T. Paul, P Leboeuf, N. Pavloff, K. Richter, and P Schlagheck, 

Phys. Rev. A 72, 063621 (2005). 
[40] M. Modugno, Phys. Rev. A 73, 013606 (2006). 
[41] N. Bilas and N. Pavloff, Eur. Phys. J. D 40, 387 (2006). 
[42] L. Sanchez-Palencia, Phys. Rev. A 74, 053625 (2006). 
[43] P. Lugan, D. Clement, P. Bouyer, A. Aspect, M. Lewenstein, 

and L. Sanchez-Palencia, Phys. Rev. Lett. 98, 170403 (2007). 
[44] G. M. Falco, A. Pelster, and R. Graham, Phys. Rev. A 75, 

063619 (2007); ibid. 76, 013624 (2007). 
[45] J. E. Lye, L. Fallani, M. Modugno, D. S. Wiersma, C. Fort, and 

M. Inguscio, Phys. Rev. Lett. 95, 070401 (2005). 
[46] D. Clement, A. F. Varon, M. Hugbart, J. A. Retter, P. Bouyer, 

L. Sanchez-Palencia, D. M. Gangardt, G. V. Shlyapnikov, and 

A. Aspect, Phys. Rev. Lett. 95, 170409 (2005). 
[47] C. Fort, L. Fallani, V. Guarrera, J. E. Lye, M. Modugno, 

D. S. Wiersma, and M. Inguscio, Phys. Rev. Lett. 95, 170410 

(2005). 

[48] J. D. Reppy, J. Low Temp. Phys. 87, 205 (1992). 

[49] G. S. Astrakharchik, J. Boronat, J. Casulleras, and S. Giorgini, 

Phys. Rev. A 66, 023603 (2002). 
[50] P Navez, A. Pelster, and R. Graham, Appl. Phys. B 86, 395 

(2007) . 



12 



[51] V. I. Yukalov and R. Graham, Phys. Rev. A 75, 023619 (2007). 

[52] V. I. Yukalov, Laser Phys. 16, 511 (2006). 

[53] N. N. Bogolubov, Lectures on Quantum Statistics (Gordon and 

Breach, New York, 1967), Vol. 1. 
[54] N. N. Bogolubov, Lectures on Quantum Statistics (Gordon and 

Breach, New York, 1970), Vol. 2. 
[55] K. Binder and A. P. Yang, Rev. Mod. Phys. 58, 801 (1986). 
[56] V. 1. Yukalov and E. R Yukalova, Phys. Part. Nucl. 31, 561 

(2000). 



[57] V. 1. Yukalov and E. R Yukalova, Phys. Part. Nucl. 35, 348 

(2004) . 

[58] V. 1. Yukalov, Phys. Rev. B 71, 184432 (2005). 
[59] S. Gluzman, V. 1. Yukalov, and D. Sornette, Phys. Rev. E 67, 
026109 (2003). 

[60] V. I. Yukalov and E. P Yukalova, Eur. Phys. J. B 55, 93 (2007). 
[61] V. I. Yukalov and E. P Yukalova, Phys. Rev. A 72, 063611 

(2005) . 



